.. _赝势说明:

赝势说明
************************************

DS-PAW目前可支持 ``.paw`` 、 ``.potcar`` 、 ``.pawpsp`` 3种格式的赝势使用，用户可通过参数 ``sys.pseudoType`` 指定使用赝势的类型。



.. _ 内部赝势:

hzw内部paw赝势
==========================================================================

DS-PAW默认使用的是 **hzw** 赝势（.paw），对应参数 ``sys.pseudoType`` 为 **-1** ，此时DS-PAW会从安装路径 **/pseudopotential** 读取赝势文件，目前 **hzw** 赝势库包含元素共72种，为元素周期表中 **1-86** 号元素（镧系元素目前只支持元素镧）。

关于 **hzw** 赝势的准确性说明：快速入门与应用案例从多个功能出发进行计算，其计算结果与文献都能很好的吻合，这验证了 **hzw** 赝势在各个功能计算中都能展现较高的准确性。

此外，针对赝势库中的 **72** 种元素所对应的单质，基于1.0版本赝势另进行了 *状态方程拟合获取稳态晶胞体积* 的计算，测试对象包括 **72** 种元素对应LDA和PBE泛函赝势，共计 **144** 个赝势文件。测试所得体积与量子化学软件 **WIEN2k** 计算所得结果进行比较，两者计算误差在元素周期表中展示如下：






.. figure:: ../imgs/phase3/pptest-LDA.png
   :align: center


.. figure:: ../imgs/phase3/pptest-PBE.png
   :align: center

| **WIEN2k 数据来源** ：https://github.com/abinit/pseudo_dojo
| **因WIEN2k网站未提供La、At计算数据，上表未展示La、At对比数据**


通过对比可得，1.0版本赝势状态方程拟合得到稳态体积与 **WIEN2k** 软件结果基本一致，其中误差最大的为元素Zn，对应的LDA和PBE赝势的误差分别为 **2.58%** 和 **3.31%** ，对这两种元素的赝势优化正在进行中。其余元素的误差基本控制在 **0.1%** 以内 ，该测试进一步验证了赝势库的整体准确性。


2024/12/31发布1.1版本赝势，稍后将公布新版赝势性能表现，敬请期待...

.. _ vasp:

VASP赝势
==========================================================================

DS-PAW提供了外部potcar格式赝势（.potcar）的使用接口，对应参数 ``sys.pseudoType`` 为 **10** ，此时DS-PAW会从默认路径 **./** 读取赝势文件。**由于版权限制，DS-PAW只提供VASP赝势的使用接口，赝势文件需用户自行准备。**
在使用时需要用户把赝势文件名做相应的修改，如果使用硅元素的LDA赝势，对应的POTCAR需要改成 Si_LDA.potcar，放置到指定目录（可以通过sys.pseudoPath设置）。

.. _ gbrv:

gbrv赝势
==========================================================================

DS-PAW提供了外部gbrv格式赝势（.pawpsp）的使用接口，对应参数 ``sys.pseudoType`` 为 **11** ，此时DS-PAW会从默认路径 **./** 读取赝势文件。gbrv赝势是一套可以免费使用的赝势库，总共提供64种元素的赝势文件，获取网站 http://www.physics.rutgers.edu/gbrv/，下载时需要注意，DS-PAW支持的格式是 PAW format for Abinit。 在使用时需要用户把赝势文件名做相应的修改，如果使用硅元素的LDA赝势，对应的赝势文件需要改成 Si_LDA.pawpsp，放置到指定目录（可以通过sys.pseudoPath设置）。




.. _ compare:

compare赝势对比
==========================================================================


.. _ 2.2.1 Si能带对比:

:math:`Si` 体系能带计算
------------------------------------------------------

为对比三种赝势的计算结果，本节以Si体系为例，分别使用以上三种赝势进行了能带计算，下图为能带对比图，该对比图说明三种赝势在描述Si的能带时，体系出较高的一致性，验证了hzw赝势的准确性。

.. figure:: ../imgs/phase3/pptest-compare.png
   :height: 600px
   :align: center

**数据来源：该图计算结果由鸿之微合作者提供**



.. _ 2.2.1 多体系带隙计算:

multi多体系带隙计算
------------------------------------------------------

本节以多个体系为例，使用1.0版本赝势对这些体系进行HSE06能带计算，并得到了多个体系的带隙值，将 **DS-PAW** 计算得到的带隙值与文献提供的 **VASP** 及 **RESCU** 软件的计算数据进行对比，得到下图，该图进一步验证了hzw赝势的准确性。

.. figure:: ../imgs/phase3/pptest-gap.png
   :height: 600px
   :align: center



**数据来源** ：https://doi.org/10.1103/PhysRevB.97.075139